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Abstract. Globular clusters house a population of compact binaries that will 
be interesting gravitational wave sources for LISA. We provide estimates for 
the numbers of sources of several categories and discuss the sensitivity of LISA 
to detecting these sources. The estimated total number of detectable sources 
ranges from ~ 10 to ~ 10 3 with gravitational wave frequencies above 1 mHz. 
These sources are typically undetectable by any other means and thus offer 
an opportunity for doing true gravitational-wave astronomy. The detection of 
these sources would provide information about both binary star evolution and 
the dynamics of globular clusters. 



1. Introduction 

Binary stars will be the most common continuous sources of gravitational radiation 
(GR) for LISA. The detection of these sources is enhanced for systems with short 
orbital periods and large component masses, favoring double compact stars, i.e., 
binaries containing white dwarf (WD), neutron star (NS), or black hole (BH) 
components. Globular clusters will be excellent hosts as they are known to 
contain large populations of compact objects in binaries formed through dynamical 
interactions (e.g., Bailyn 1995; Hut et al 1992). While some of these systems have 
been detected directly through their electromagnetic radiation (e.g., as radio pulsars, 
X-ray or UV sources), most are thought to be invisible in conventional astronomy, 
but they could be very bright in gravitational radiation. For example, a large, dense 
globular cluster such as 47 Tuc may contain as many as ~ lfj 4 NS, of which ~ lfj 3 
may be active radio millisecond pulsars, but only ~ 20 are bright enough to have been 
detected in pulsar surveys (Camilo et al 2000). Most of the detected pulsars are in 
binaries with low-mass WD companions. 

If LISA can detect and identify some of these sources within specific globular 
clusters, then gravitational wave astronomy will be able to contribute significantly 
to our understanding of binary star evolution and stellar dynamics in dense cluster 
systems. Many Galactic globular clusters are close enough that GR emission from 
double degenerate binaries may be detected by LISA with a signal-to-noise (S/N) 
ratio > 2. At this level, LISA will have an angular resolution of a few square degrees. 
Individual clusters will be smaller than this, but the typical separation between clusters 
is much larger, so that LISA will be able to separate sources in different clusters. 
One of the primary sources of noise in the frequency band of cluster binaries is the 
Galactic disk population of close white-dwarf binaries, which is expected to produce a 
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confusion- limited signal above the noise level of LISA (Hils and Bender 2000, Schneider 
et al 2000). It may be possible to avoid this noise by looking specifically for signals 
originating out of the plane of the Galaxy, such as those arising from globular clusters. 
In this paper, we discuss the possible source populations and we present some simple 
methods for crudely estimating their numbers. 

Since the angular resolution of LISA depends upon measuring the Doppler shift 
in frequency as a phase variation during the motion of the detector about the sun, 
the wavelength of the signal must be less than the size of LISA's orbit. This places 
an upper limit on the binary orbital period, P or b = 2// gw = 2A/c < AAU/c ~ 2000 s. 
In addition, we require that the S/N ratio be > 2 for the binary to be detectable 
with the required positional accuracy. By coincidence, this also imposes a maximum 
period P max — 2000 s above which the confusion-limited noise from Galactic binaries 
rises significantly. Finally, we require that the systems be detached so that their period 
evolution is solely determined by GR, providing a "clean" signal. This should also 
improve the extraction of the signal from the confusion-limited background. These 
requirements limit the viable sources to double compact objects. We will not discuss 
other standard LISA sources such as contact binaries containing two main-sequence 
stars or a main-sequence star and a WD, as their signals will not be strong enough 
or clean enough. Integrating the in-spiral of a binary driven by GR alone from P max 
down to contact, we can write the total length of time spent in the LISA band as 

/ M u- \" 5/3 /P \ 8/3 
Tlisa — 4 x 10 6 yr ( _2£i±l£E ) ( £™L ) (1) 

where M c hi rp = (mim 2 ) 3, ' 5 /(' n i + m^) 1 ^ is the chirp mass, given in terms of the 
binary component masses m\ and m^- 

2. Theoretical Estimates 

Numerical simulations of globular cluster dynamics can be used to try to estimate 
the numbers and properties of different classes of binary systems from first principles. 
Here we illustrate this approach from three recent studies using different numerical 
methods. 

2.1. N-body simulations 

The dynamical simulations of Portegies Zwart and McMillan (2000, hereafter PZM) 
were developed to estimate BH-BH merger rates for LIGO and VIRGO. In their TV- 
body simulations of globular cluster evolution, a component of 10 M© primordial BH 
quickly undergoes mass segregation, forming an effectively isolated subsystem at the 
cluster core in ~ 10 9 yr. Through dynamical interactions, hard BH binaries are then 
formed and ejected from the cluster and the entire BH subsystem evaporates typically 
within ~5x 10 9 yr. The basic conclusion of this work is that globular clusters today 
should no longer contain any of their primordial BH. This explains naturally why no 
BH X-ray binary has ever been detected in a globular cluster, even though clusters 
contain many bright NS X-ray binaries (see, e.g., Deutsch et al 2000). A large number 
of ejected BH binaries will coalesce within a Hubble time, and these may provide a 
significant event rate for ground-based interferometers. 

However, the last BH binary remaining in the cluster core may not be ejected, 
since it can now only interact with the (much less massive) normal stars. The 
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dynamical simulations of PZM suggest that the initial separations of typical BH 
binaries should be <~ 30 — 300 Rq. The eccentricities are typically large, and roughly 
half of these binaries will coalesce through GR within the remaining t coa i ~5x 10 9 yr 
of the cluster evolution to the present (it is easy to show that, for coalescing 
binaries, further hardening through dynamical interactions with normal cluster stars 
is negligible, and that ejection from the cluster is extremely unlikely). Unfortunately, 
with Tlisa — 10 5 yr for two 10 Mq BH (see eq. 1), even if we assume optimistically that 
~ 100 globular clusters in our Galaxy retained a coalescing BH binary, we still conclude 
that the number of BH-BH systems in the LISA band today is <~ 100 TusAAcoai <C 1. 
Thus BH-BH binaries in clusters are not promising sources for LISA. 



2.2. Encounter Rates 

Davies (1995) calculated the production rates of various binaries through dynamical 
interactions using a combination of scattering experiments and static models of 
globular clusters. His results can be used to estimate the numbers of NS-NS, NS- 
WD, and WD- WD binaries with P orb < 2000s (Benacquista 1999). 

The original dynamical simulation calculates three-body interaction cross-sections 
for a variety of interactions and outcomes and then determines encounter rates for a 
number of model globular clusters. Production rates of binaries were inferred from the 
cross-sections applied to the encounters. As a side result, the merger rates through 
GR induced in-spiral of binaries were produced for 15 Gyr runs. The runs were also 
extended until all binaries had merged. This generally took another 15 Gyr. With 
these two results, a merger rate can be determined through 

V= 15 Gyr (2) 

where is the number of mergers after all binaries have merged and 7115 is the 
number which have merged after the first 15 Gyr. The evolution of the orbital period 
can be determined from P = — fc Q P~ 5 / 3 with 

O C 

Combining these results in a number distribution given by 

dn = ^-P 5 / 3 dP. (4) 

Using the numbers found in Davies (1995) and scaling up to the entire globular cluster 
system, with a total number of stars N gc ~ 10 7 5 , we find 16 NS-NS binaries, 125 NS- 
WD binaries, and 176 WD- WD binaries in the LISA band. 



2.3. Monte Carlo Simulations 

Dynamical Monte Carlo simulations can be used to study self-consistently the 
evolution of binary populations within evolving globular cluster models (Rasio 2000). 
Rasio et al 2000 used this approach to study the formation and evolution of NS- 
WD binaries, motivated by the recent detection of a large number of millisecond 
radio pulsars with WD companions in 47 Tuc (Camilo et al 1999). The simulations 
include an evolving core of the background cluster, stellar evolution, mass segregation, 
and dynamical 3-body interactions. In the process of generating the appropriate 
populations of binary millisecond pulsars, the simulations also produce populations of 
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NS-WD binaries in the period range of interest for LISA. There are two such binaries 
in the simulation of 47 Tuc (Rasio et al 2000) . From this admittedly small sample it 
is possible to place an order of magnitude estimate on the number of such binaries 
in the globular clusters system. Assuming that the 2 NS-WD binaries (presently in 
the LISA band) for the ~ 10 6 M Q in 47 Tuc is consistent throughout the globular 
cluster system, we would predict a total of ~ 60 such binaries for the ~ 10 7 ' 5 M Q in 
all clusters. Note that this number is within a factor of two of the number estimated 
above from encounter rates. 

3. Semi-Empirical Estimates 

Alternatively we can estimate the numbers of the various possible LISA sources by 
trying to relate them to certain types of binaries already observed in electromagnetic 
radiation. Here the best candidates are NS-WD systems, also observable as X-ray 
binaries, and NS-NS systems, also observable as binary radio pulsars. 

3.1. NS-WD Binaries 

The small number of observed ultra-compact X-ray binaries (UXBs) in globular 
clusters (Dcutsch et al 2000) can be used to estimate empirically the number of 
NS-WD binaries in the LISA band. UXBs are bright enough that we essentially see 
all of them in the Galactic globular cluster system. Assuming that they originate from 
the GR-driven in-spiral of NS-WD binaries all the way to contact (Rasio et al 2000) 
and that the process is ongoing, we can estimate the number of progenitor systems in 
the LISA band by looking at the time spent by a system in both phases. Let Ax be 
the number of UXBs, and Tx their typical lifetime. Also, let -/Vlisa be the number of 
progenitors in the LISA band, with the corresponding Tlisa given by eq. 1. Then, if 
the process is stationary, we must have 

N x A L i S a 

7^ 7^ • (5) 

JX J- LISA 

Of the 12 known bright X-ray binaries in clusters, 3 are known to have very short 
orbital periods (< 1 hr) and 3 are known to have long orbital periods (Deutsch et al 
1999). Thus, 3 < Nx < 9. The lifetime Tx is rather uncertain, depending on a number 
of poorly understood processes. A standard treatment of mass transfer driven by GR 
only would give Tx ~ 10 9 yr, but other effects such as irradiation or tidal heating of 
the companion may well shorten this to Tx ~ 10 7 yr, which would then also provide 
consistency with the birthrate of binary millisecond pulsars (Applegate and Shaham 
1994; Tavani 1991; Rasio et al 2000). For a typical system containing a 1.4 M NS and 
a 0.5 Mq WD companion, eq. 1 gives Tlisa — 7 x 10 6 yr. Adopting the most optimistic 
values of Ax and Tx we arrive at Alisa ~ 1 — 10, about an order of magnitude smaller 
than predicted by dynamical simulations. 

3.2. NS-NS Binaries 

One NS-NS binary (M15C) has been detected directly as a binary radio pulsar, out of a 
total of ~ 40 radio pulsars known in all Galactic globular clusters. A proper statistical 
analysis taking into account the various selection effects in pulsar surveys is well 
beyond the scope of this paper. However, since the dynamical processes that created 
M15C are thought to happen quite generically in all dense clusters (see, e.g., Phinney 
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and Sigurdsson 1991), and since selection effects, if anything, tend to make pulsars in 
close binaries more difficult to detect than other types of radio pulsars (Johnston and 
Kulkarni 1991), it may not be unreasonable to estimate the total number of NS-NS 
binaries in globular clusters as ~ JVns/40, where Ans is the total number of NS in 
clusters. A strict lower limit on A%s comes from estimates of the total number of active 
radio pulsars in clusters, giving A^ns > 10 4 (see, e.g., Kulkarni et al 1990). An upper 
limit comes from dynamical models of globular clusters suggesting that as much as 
~ 5% of the total mass of a cluster may be in the form of 1.4 M© dark remnants (some 
of which may be massive WDs, which, for our purposes, are indistinguishable from 
NS; see, e.g., Murphy et al 1998). Taking again the total mass in globular clusters 
~ 10 7 ' 5 M© (Prince et al 1991), this gives N?is < 10 6 . Adopting the coalescence time 
of M15C, i CO ai = 2 x 10 8 yr, as typical, and using Tlisa = 3 x 10 6 yr from cq. 1 for two 
1.4 M© NS, we get a number ~ 1 - 100 NS-NS binaries in the LISA band, in rough 
agreement with the predictions of dynamical simulations. 



4. Detection Efficiency 

From the above estimates, we see that there may be many compact binaries in the 
Galactic globular cluster system that will be detectable GR sources for LISA. In 
order to determine the efficiency of LISA for detecting these sources throughout the 
globular cluster system, we calculated the S/N ratio for a range of binary masses and 
frequencies from each globular cluster in the system. We begin by modeling the strain 
amplitude measured by LISA. A stationary detector will measure a strain amplitude 
found using the point-mass quadrupolc formula for a circularized binary given by 

h(t) = ^-{A + F+ cos (2wft) + A X F X sin (2nft)), (6) 

where F + and F x represent the sensitivity of the detector to each polarization. The 
amplitudes are given by 

A + = ^M*£ p f*/*{l + a x>i) (7) 

A x =-2^M c ^ p / 2 /3 cosi . (8) 

The angle of inclination, i, is the angle that the line of sight makes with the orbital 
plane of the binary. Since LISA is in orbit about the sun and its orientation changes, 
the actual signal received by LISA is given by 

nt 

(9) 



/q 

h(t) = ^-A(t) cos 



J 2nf{t')dt' + ip p {t)+<p D (t) + <p 



The Dopplcr shift in frequency due to LISA's motion about the sun is incorporated in 
(po(t), and the changing orientation of LISA with respect to the source is incorporated 
in both tp p (t) and A{t) = ^(A+F+) 2 + (A X F X ) 2 . The exact form of <p D , ip p , F+, 
and F x can be found in Cutler (1998). The S/N ratio p is then found from 

~ h * M)ha{f) swr (10) 

with the sum over a representing the sum over both semi-independent arms of LISA 
and the Fourier transform h defined as 

/>oo 

M/)= / e 27Tift h a (t)dt. (11) 
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Using globular cluster position and distance data from Harris (1996), we 
determined the S/N ratio for NS-WD binaries with WD masses in the range of 0.1 
to 0.4 M Q and frequencies ranging from 1 mHz up to Roche lobe overflow, using the 
expression of Rappaport et al (1987), / max = 43.3 mHz (Mwd/Mq). We note that 
the signal strength for a binary with i = is about twice that for one with i = tt/2. 
The globular clusters were then classified by the strength of the detected signal for 
binaries with favorable or unfavorable angles of inclination. There are 7 clusters with 
some binaries having i — tt/2 and p > 5 and 17 for which only binaries with i = 
have p > 5. The binaries associated with these clusters would have sufficient angular 
resolution to clearly identify the binary with its cluster. There are 72 clusters with 
i = n/2 and p > 2 and 36 with only i = and p > 2. These clusters could house 
binaries whose angular resolution would be a few square degrees, and the identification 
of the binary with its host cluster would be a little more suspect. There are only 15 
clusters which are too distant to detect any binaries with p>2. 

5. Conclusion 

There is ample, direct and indirect, evidence for the existence of a population of 
compact binaries in globular clusters that will be detectable GR sources for LISA. The 
observed populations of X-ray binaries and millisecond radio pulsars clearly suggest a 
rich and interesting variety of dynamical processes that should produce large numbers 
of compact binaries. A variety of estimation techniques, both theoretical and semi- 
empirical, are in rough agreement to predict as many as ~ 100 NS-NS, NS-WD, 
and WD- WD binaries detectable by LISA in the Galactic globular cluster system. 
Although there are very large uncertainties in the exact numbers, the existence of these 
populations is not in doubt. LISA will be able to probe about 90% of Galactic globular 
clusters to search for binaries with orbital periods below 2000 s. Measurements of this 
population of binaries will provide valuable insight into the dynamical evolution of 
dense star clusters and the role that binaries play in them. 
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